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We apply multifractal analysis to an experimentally obtained quasi-two-dimensional crys- 
tal with fourfold symmetry in order to characterize the sidebranch structure of a dendritic 
pattern. In our analysis, the stem of the dendritic pattern is regarded as a one-dimensional 
support on which a measure is defined and the measure is identified with the area, perimeter 
length, and growth rate distributions. It is found that these distributions have multifractality 
and the results for the area and perimeter length distributions, in the competitive growth 
, regime of sidebranches, are phenomenologically understood as a simple partitioning process. 
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1. Introduction 

> 

(N 

It has been well established that a dendritic pattern is typical in crystal growth. Its growth 
process is dominated by diffusion and anisotropy 1 and consists of several stages; (i) A tip 
grows stably and steadily, and a straight stem is formed. The stability of the tip is attributed 
to anisotropy. (ii) Sidebranches are generated behind the tip due to noise effects 2 and the 
instability of a flat interface. 3 (hi) Sidebranches grow competing mutually. This competition 
is known as one of the most characteristic and interesting properties in growth dominated by 
diffusion. A longer sidebranch screens off the diffusional field and suppresses the growth of 
shorter ones around it. This process occurs on various length scales. As a result, complicated 
and hierarchical structures are formed, (iv) Finally, surviving branches grow independently. 

Our aim in this paper is to characterize the sidebranch structure (at stages (ii) and (Hi)) of 
a dendritic pattern. It is interesting and worth considering since the sidebranch structure 
determines the outline of the pattern. Many types of scaling analysis based on scaling idea have 
been attempted. One concerns the properties of global structure constructed by sidebranches. 
For the pattern of a three-dimensional dendritic crystal projected onto a two-dimensional 
plane, it has been reported 4-6 that S(X) ~ X &s and L(X) ~ X &L with 5 S ~ 5 L ~ 1.7, 
where S is the area of the pattern and L the perimeter length up to X, which is the distance 
along the stem from the tip. Interestingly, for the pattern of a quasi-two-dimensional crystal, 
S(X) ~ X Ss with 5,s ~ 1.5. 7 Active sidebranches, whose growth is not suppressed by the 
screening effect of a longer sidebranch, form the envelope Z{X) of the pattern, where Z(X) 
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denotes the height from the stem at X. The envelope is constructed by connecting the tips 
of active sidebranches and obeys the power law Z(X) ~ X K . 5 > 6 > 8 It has been observed that 
k < 1 near behind the tip and k > 1 far from the tip. 9 Another scaling analysis concerns the 
growth of an individual sidebranch or the statistical properties of a set of sidebranches. In the 
competing growth of sidebranches, each branch grows as hj(tj) ~ ty, with Oj ~ 0.5 - 0.7, 
where hj(tj) is the height of the j'-th branch at tj, the time from its birth, and subsequently the 
growth decays exponentially. 10 For a set of sidebranches it is found that the height distribution 
N(h) obeys the following power law, N(h) ~ h" 13 , with j3 ~ 2. 2. 11 

In this paper we present a new characterization of the sidebranch structure of a dendritic 
pattern using multifractal formalism. Since the stem of a dendritic pattern grows straight, it 
can be considered as a one-dimensional support on which a probability measure is defined. We 
apply this approach to an experimentally obtained quasi-two-dimensional dendritic pattern 
of an NH4CI crystal. The crystal is obtained from a supersaturated solution and has fourfold 
symmetry. We identify the probability measure with the area and perimeter length distri- 
butions of the pattern and the growth rate distribution at the interface. In the competitive 
growth of sidebranches, the solute particles diffusing in the solvent are distributed unequally 
to branches since it is easy to reach the tip of longer branches whereas it is difficult to reach 
the shorter branches between longer ones. This process is expected to occur on various length 
scales, similar to energy dissipation in turbulence, which is simply modeled by the "binomial 
branching process" 12 and shows multifractality. As for the growth rate distribution on the 
interface, it is known to have multifractality with the interface itself as a (fractal) support. 13 
Therefore we expect multifractality in our point of view here. 

The rest of the paper is organized as follows: In section 2 our crystallization experiment is 
briefly described. The multifractal formulation is given in section 3. The binomial branching 
process is introduced and the results of multifractal analysis are referred to. In section 4 our 
results and discussion are given. For the area and perimeter length distributions, comparison 
with those for the binomial branching process are presented. Section 5 is dedicated to the 
summary and future outlook. 

2. Experiment 

We analyze a quasi-two-dimensional dendritic pattern, with well-developed sidebranches 
and a clear envelope, obtained from an NH4CI solution growth experiment. The details of the 
experiment are described in our previous articles: 11 An NH4CI aqueous solution saturated at 
approximately 40 °C is sealed in a Hele-Shaw cell, which has a small gap between two glass 
plates placed in parallel. The thickness of the gap is 100 fim. Then when the temperature is 
lowered to approximately 30 °C, the solution becomes supersaturated and nucleation takes 
place. The direction of the tip growth is (100) in the supersaturated solution. Sidebranches 
grow perpendicularly to the stem, with small sub-sidebranches perpendicular to them. The 
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Fig. 1. Image of a dendritic crystal. The resolution is 640x480 pixels. For our analysis, the X-axis is 
set along the stem with the tip on it. The Z-axis is set parallel to the sidebranches. The pattern 
is covered with strips of width e aligned parallel to the Z-axis. There is no correlation between 
the sidebranch patterns of Z > and Z < 0. 



observed tip velocity va p is 40—49 ^m/sec. The diffusion length of the tip, Zd = 2D/w t i p , 
where D is the diffusion constant of NH4CI (2.6 x 10 3 ^m 2 /s 14 ), is larger than the thickness 
of the cell. Therefore the growth is considered to be quasi-two-dimensional. The image of the 
crystal is obtained by using a microscope and charge-coupled device (CCD) camera and is 
binarized by an image processing. The image of a crystal is shown in Fig. 1, whose resolution 
is 640 x 480 pixels. 

3. Multifractal 

3.1 Formulation: on one- dimensional support 

Suppose that the stem of a dendrite grows along the X-axis, with the tip on the X-axis, 
and the sidebranches grow in the positive and negative Z directions (see Fig. 1). There is 
no correlation between the growths in the positive and negative Z regions. 2 Therefore each 
pattern can be dealt with as an independent sample. 

Consider that a pattern (for example, the pattern shown in Fig. 1, Z > 0) is covered with 
disjoint strips of width e aligned parallel to the Z shown in FIG.l. Let Pj(e) be a 

measure (nonnegative scalar quantity) assigned to the j-th strip. The measure is normalized 
to be a probability 

JV(6) 

X>(*) = 1, 

3=1 
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where N(e) is the number of strips necessary to cover the pattern completely. In our point of 
view, the stem is regarded as a one-dimensional support on which the probability measure is 
defined. 

The partition function Z(q, e) is defined by the probability measure as 

Z(q,e) = £ [ P M q - (2) 

Based on the expectation that for small e the scaling law Z(q, e) ~ e T ^ holds, the multifractal 
exponent r{q) is defined as 

T (q) = ]im — . (3) 

e->o log e 

Practically r{q) is evaluated from the slope of log Z(q, e) versus log e. Then the generalized 
dimension D(q) is defined as 15 

D(q) = ^jr(q), for q±\, (4) 
D(l) = 

= ^S^jggg^, for q = 1 . (5) 
e^o loge 

Using the Legendre transformation the singularity exponent a and its fractal dimension f(a) 
are given as functions of q as 16 

«( 4 ) - ^ <«> 

/(«(?)) = q<*(q)-T(q). (7) 

However, it is not useful to numerically evaluate a and f(a) from Eqs. (6) and (7), since they 
may produce relatively large errors. Therefore instead, we adopt the direct method described 
below. 17 

Let us construct a new probability measure /j,j(e,q) with parameter q from Pj(e) as 

*(«■«) = J^ } ' . w m 

E j {Pj[£)} q 

Then let us define £(e, q) and £(e, g) as 

CM) = 2/*j(e,?)Iogbj(e)]. (9) 

= $^Atj(e,g)Iog[/ij(e,g)]. (10) 
i 

From them a and f(a) are given as functions of q as 

a(g) = lim^^, (11) 
e->o loge 
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f{q) = lim^^. (12) 
JyHJ e^o loge K J 

Practically they are evaluated from the slopes of ((e,q) and £(e, q) versus loge, respectively. 
Direct calculation shows that the definitions (11) and (12) satisfy the relations (6) and (7). 



3.2 Binomial branching process 

We refer to the binomial branching process 12 for our later consideration. It shows mul- 
tifractality and fortunately the multifractal spectrum can be exactly calculated due to its 
simplicity. 

Suppose that a segment of length 1 is divided into two segments of length 1/2. A prob- 
ability measure p (> 1/2) is assigned to the left segment and (1 — p) to the right. This p 
is the only adjustable parameter of the process. Next, each segment is subdivided into two 
equal halves and the measure is partitioned into p to the left half and (1 — p) to the right. 
So there are four segments of length 1/4 and the measures p 2 , p{l —p), (1 — p)p and (1 — p) 2 
are assigned to the segments from left to right. This procedure is repeated again and again 
(see Fig. 2, we can see a similarity between the n = 8 pattern and the sidebranch structure 
shown in Fig. 1). At the n-th iteration, there are 2™ segments of length 2~ n and the number 
of segments with measure p k (l —p) n ~ k , k = 0, 1, • • • , n, is Qf) = n\/[k\(n — k)\]. Therefore the 
partition function of the stage, Z(q,n) = Z(q,e) where e = 2 _n , is immediately calculated as 

z{ q ,n) = e( ? Va-prY 

k=0 \ k / 

= |p9 + (l-p)9]». (13) 

Then the exponent r(q) can be obtained as 

w+v-M (14) 

log 2 

Using the Legendre transformation the singularity exponent and the fractal dimension are 
calculated as 

a(q) = _V^P+(l-ri)log(l-p) 
y ' log 2 y ' 

f {<*(«)) = - r/l0gr? + (1 l0 "^ )l0g(1 " ??) , (16) 

where n = p q /[p q + (1 — p) q \. The direct evaluation using eqs.(ll) and (12) gives the same 
result. 

The f(a) spectrum takes a continuous value for [a m i n , a max ], where a m j n = — log 2 p and 
"max = — log 2 (l — p). It is symmetric with respect to a = a(q = 0) = — (log 2 p + log 2 (l — p))/2 
and takes the maximum at a = a(q = 0), f(a(q = 0)) = D(0) = 1, reflecting the fact that 
the support is one-dimensional. The information dimension D(l) = (a(q = 1) = f(a(q = 1))) 
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Fig. 2. (a) Different stages of the binomial branching process. Each segment is divided into two equal 
subsegments at the next stage. The measure is divided into nonequal fractions, p and (1 — p). (b) 
f(a) spectrum of the binomial branching process with p = 0.7. 



is given as 



D(l) 



plogp + (1 -p) log(l -p) 
log 2 



(17) 
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Fig. 3. Generalized dimension D{q) of the area distribution for the pattern shown in Fig. 1, Z > 0, 
and that for the binomial branching process with p = 0.68. The error bars are obtained from the 
data of 30 samples. 



4. Results and Discussion 

4-1 Area distribution 

First of all, we investigate the structure of the area distribution. The area of the pattern 
is defined as the number of pixels which constitute the pattern. 

The generalized dimension D(q) for the pattern shown in Fig. 1, Z > is shown in Fig. 
3. Taking into consideration that each sidebranch has a finite thickness of ~5— 25 pixels, we 
use 6 pixels as the minimum of the strip width in our analysis. On the other hand we use 120 
pixels as the maximum. The exponent r(q) is obtained by least squares method. The log-log 
plots of Z(q, e) against the strip width e and the fitting lines for some values of q are shown 
in Fig.4. The scaling relation holds quite well. 

The multifractal f-a spectrum for the same pattern is shown in Fig. 5, along with that 
for the binomial branching process. The maximum of f(a) equals 1, since the support is one- 
dimensional. The spectrum is evaluated using eqs. (8) — (12). The plots of C(<L e ) an d e) 
against log e and least squares fitting are shown in Fig. 6. The fitting gives considerably good 
evaluations, although there exist some scattering points for larger q. 

It is found that both D(q) and the f-a spectrum take continuous values dependent on q 
and a, respectively, within the range between D(q — >■ oo) = a m \ n ~ 0.56 and D(q — v — oo) = 
"max ~ 1-7. This finding indicates that the area distribution has multifractality. Multifractality 
characterizes the pattern in more detail than the global scaling relations mentioned in the 
Introduction by using a local scaling relation and distributions. That is, under the condition 
of the global scaling relation S(X) ~ X Ss (see the Introduction), at around a certain point 
the area scales locally as ~ e a (e — >■ 0) and such a point is distributed as ~ e _ ^( Q ) . We compare 
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Fig. 4. (Color online) Log-log plots of Z(q, e) of the area distribution for the pattern shown in Fig. 
1, Z > 0, against strip width e, for q —0, 1, 2, and 3. The slope gives the exponent r(q). For the 
ease of viewing, the measure is not normalized to be a probability, i.e., Z(q = 1, e) is the total 
area of the pattern. 




the results with the binomial branching process with p = 0.68, which gives the same a m i n . 
There is a good agreement in the positive q or small a region. 

The result is almost independent of the pattern size, as long as the pattern is large enough 
that competitive growth of sidebranches is well-developed. To see this, in Fig. 7 we show the 
f(ot) spectra for the full pattern shown in Fig. 1, Z > 0, and for the partial patterns from 
the tip to 75%, 66%, 33%, and 25% of the full pattern. There is little difference between the 
spectra of the three partial patterns of 100%, 75%, and 66%. These can be considered as 
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Fig. 6. (Color online) Fitting of a and / for q =0, 1, 2, and 3. (a) Plots of £(<7, e) of the area 
distribution for the pattern shown in Fig.l, Z > 0, and (b) plots of £((?,e) against lne. 



sufficiently large samples. However, for the 33% and 25% patterns, the minimum singularity 
exponent a m ; n is larger than those for the former three patterns. Therefore the small-a part 
of the spectrum is contributed to by well-developed longer sidebranches and the large-a part 
by small sidebranches near the tip or deep inside the forest of longer sidebranches. 

Our results of some characteristic exponents are summarized in Table. I, along with those 
for the binomial branching process with p = 0.68. They show a good agreement, except for 
c^max- Therefore, our phenomenological scenario mentioned in the introduction seems reason- 
able in the competitive growth regime. However, how the unequal distribution of solute to 
sidebranches on various length scales is related to the diffusional growth mechanism and how 
the value p = 0.68 is derived still remain unclear. On the other hand, there exists some dis- 
agreement slightly to the right of the peak of the spectrum This is attributed to the fact that 
near behind the tip sidebranches are short and growing almost independently, not competing 
with each other. Note that the reliability of the spectrum for q < or large a is considerably 
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Fig. 7. (Color online) Size dependence of the f-a spectrum of the area distribution for the patterns 
shown in Fig.l, Z > 0. 



"min a max a(q = 0) D(l) 

Area 0.54±0.06 1.58±0.11 1.05±0.02 0.94±0.03 
p=0.68 0.56 1.64 1.10 0.90 

Table I. List of characteristic scaling exponents for the area distribution, along with those for the 
binomial branching process with p = 0.68. Note that a(q = 0) is the value at which the fractal 
dimension takes maximum / = 1. Data are obtained over 30 samples. 



lower than that for q > or small a due to the limitation of the resolution. 

4-2 Perimeter length 

The perimeter length is defined as the number of pixels which constitute the interface of 
the pattern. A pixel is said to constitute the interface if it is a part of the pattern and at least 
one of its four neighboring pixels is not. 

Figure 8 shows the generalized dimension D{q) and the f-a spectrum for the pattern 
shown in Fig.l, Z > 0. It is clear that the distribution has multifractality. Figures 9 and 10 
show the log-log plots of Z(q,e) against e and the plots of ((q,e) and e) against loge, 
respectively. They show that the scaling relation holds well. 

Some characteristic exponents are summarized in Table. II. Note that for some samples 
/(a max ) does not equal zero, as shown in Fig. 8. For the perimeter length distribution, the 
number of pixels used in the calculation is considerably smaller than that for the area distri- 
bution, since only the pixels constituting the interface of the pattern are to be studied. Due 
to this fact and the limitation of the resolution, it is quite difficult to obtain a result with 
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Fig. 8. (a) Generalized dimension and (b)/-a spectrum of the perimeter length distribution for the 
pattern shown in Fig. 1, Z > 0. The error bars are obtained from the data of 30 samples. 




q=0 
q=l 
q=2 
q=3 



Fig. 9. (Color online) Log-log plots of Z(q, e) of the perimeter length distribution for the pattern 
shown in Fig. 1, Z > 0, against the strip width e, for q =0, 1, 2, and 3. The measure is not 
normalized, i.e., Z(q = 1, e) is the total perimeter length of the pattern. 



11/17 



J. Phys. Soc. Jpn. 



Full Paper 



-l 
-1.5 
-2 
-2.5 

dr 

& -3.5 
-4 
-4.5 
-5 
-5.5 

1 1.5 2 2.5 3 3.5 4 4.5 5 

-1 
-1.5 
-2 
-2.5 
S -3 
-3.5 
-4 
-4.5 
-5 

1 1.5 2 2.5 3 3.5 4 4.5 5 

Fig. 10. (Color online) Fitting of a and /. (a) Plots of ((q, e) of the perimeter length for the pattern 
shown in Fig. 1, Z > 0, and (b) plots of £(q, e) against lne for q —0, 1, 2, and 3. 

satisfactory precision for the q < or large-a region and to conclude whether /(a max ) ^ 
or not. It is found that a m i n for the perimeter length distribution is larger than that for the 
area distribution. However, there is a good agreement for the values of a(q = 0) and D{1) 
with relatively small error. This is interpreted as follows: a m ; n is dominated by the contri- 
bution from long and thick sidebranches only, while both from thick and long, and thin and 
short sidebranches contribute to a(q = 0) and D(l). For thick sidebranches the difference 
between the number of pixels constituting the branch and that constituting the interface is 
large, but for thin branches this difference is small. Therefore it is reasonable to assume the 
same scenario corresponding to the binomial branching process for thin branches as for the 
area distribution. 

4-3 Growth rate distribution 

In principle, it is a faithful method to the original data to evaluate the growth rate from 
the growth area between two successive images. However, for a dendritic pattern such a 
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"min "max a(q = 0) D(l) 

Perimeter 0.67±0.07 1.28±0.16 1.05±0.01 0.95±0.01 

Table II. List of characteristic scaling exponents for the perimeter length distribution. Data are 
obtained over 30 samples. 



method is quite difficult to implement with satisfactory precision due to the limitation of the 
resolution and since the difference of the growth rates between in the fast region and in the 
slow region is quite large. Therefore instead, we evaluate the growth rate p(r) at point r on the 
interface by numerically solving the Laplace equation V 2 ^>(r) = outside of the pattern, where 
0(r) denotes the concentration field, on a square lattice. We set, as a boundary condition, 
0(r) = Const., uniformly on the interface and evaluate the growth rate as the gradient of the 
concentration field: 

p(r) ~ |V0(r)|. (18) 

This evaluation is valid if the diffusion length is larger than the characteristic length of the 
system - the tip radius of the stem or the average spacing of sidebranch generation - and 
this is the case: the diffusion length is longer than 100 fim and the characteristic length is of 
the order of 1 fim. We neglect the surface tension effect, since the surface tension effect does 
not affect the result in the unscreened large growth region. 13 We are intersted in the scaling 
structure in that region. Then the measure Pj(e) in the j-th strip is given as 

rgj— th strip 

which is normalized to be a probability. 

Figure 11 shows the generalized dimension for q > and the multifractal f-a spectrum 
for small a of the pattern shown in Fig. 1, Z > 0. The log-log plot of Z(q,e) against e 
and the plots of ((q,e) and e) against loge are shown in Figs. 12 and 13, respectively. 
Multifractality of the distribution and scaling property are observed. Some characteristic 
exponents are summarized in Table. III. It is unclear how these values are derived. However, 
the growth rate distribution considered here is the harmonic measure redefined on a one- 
dimensional support, which is originally defined on a fractal interface embedded in the two- 
dimensional plane. Therefore it is expected that our results are closely related to the results 
for the harmonic measure given on the interface of a dendritic pattern. 13 Particularly it is 
known that the information dimension for the latter, Df(l), is exactly proved to be l. 18 We 
conjecture that there will be a simple relation between the two values of the information 
dimension, D(l) = Df(l)/Df(0) = 1/Df(0) where D^(0)(~ 1.5) is the fractal dimension of 
the dendrite interface. This conjecture is based on the assumption that since we consider the 
same measure on different supports, the interface with the fractal dimension Df(0) and the 
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Q!min a(q = 0) D(l) 
0.3±0.1 1.6±0.2 0.70±0.06 

Tabic III. List of characteristic singular exponents for the growth rate distribution. Data are obtained 
over 30 samples. 



stem with D(0) = 1, the two multifractal spectra corresponding to these supports are similar 
and the condition that they contact the line f(a) = a. Our result is consistent with the 
conjecture within error. 

The small-a region is contributed to from unscreened active growth. In our situation 
there are two domains where growth is active: one is around the tip of longer sidebranches 
and the other is near behind the tip of the stem, where small short branches are growing 
almost independently. Therefore, the results for the growth rate distribution are different from 
those for the area and perimeter length distributions due to the difference of the structure of 
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12. (Color online) Log-log plots of Z(q, e) of the growth rate distribution for the pattern shown 
in Fig. 1, Z > 0, against the strip width e, for q =0, 1, 2, and 3. The measure is not normalized, 
i.e., Z(q = 1, e) is the total sum of the growth rates. 
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contribution. For example, short branches contribute to the small-a part of the spectra of the 
area and perimeter length distributions. On the other hand, for the growth rate distribution, 
active short branches near behind the tip contribute to the small-a part but frozen ones deep 
inside the forest of longer branches to the large-a part, which is not discussed here. 

5. Summary and Outlook 

In order to characterize the sidebranch structure of a dendritic pattern we applied multi- 
fractal formalism to the fourfold pattern of an NH4CI quasi-two-dimensional crystal, regarding 
the stem as a one-dimensional support on which the probability measure is given. Multifrac- 
tality of the area and perimeter length distributions was manifested and was well understood 
phenomenologically as the binomial branching process. Furthermore the growth rate distri- 
bution also showed multifractality from the point of view on one-dimensional support. 

There are some problems for further understanding. One concerns the relation between 
diffusional the growth mechanism and the binomial branching process, which we use to phe- 
nomenologically understand the result. It will provide a strong support to our consideration 
if it is clarified. Another concerns higher multifractality. Since the growth rate is originally 
given on the interface, it seems natural to consider a measure which has multifractality and 
is defined on a support which also has multifractality. To our knowledge, the study of such 
a system has not been carried out in detail, except for the cases of certain simple models. 19 
In this point of view the stem is considered as a one-dimensional base space on which the 
support is distributed. At a certain point on this one-dimensional base space, it is expected 
that the support locally scales as ~ e ai and the measure also locally scales as ~ e a ' 2 . Such 
a point is distributed as ~ e -- F ( Q i. Q 2) ) w here F{a\,a2) is a function of both a\ and a<i- We 
hope that useful information will be derived from F(a\, 0:2), which will enable a more detailed 
understanding of the sidebranch structure. 
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